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Abstract 

Directed flow of antiprotons is studied in Au+Au collisions at a beam momentum of 
11. 5A GeV/c. It is shown that antiproton directed flow is anti-correlated to proton 
flow. The measured transverse momentum dependence of the antiproton flow is 
compared with predictions of the RQMD event generator. 
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1 Introduction 

The production of antiprotons at AGS energies (10-15A GeV) is near thresh- 
old in nucleon-nucleon coUisions. Via collective effects, p production can be 
enhanced in nucleus-nucleus collisions as compared to p-|-p collisions [1-4]. 
The observed yield of antiprotons is a result of both production and subse- 
quent annihilation. The high antiquark densities in the quark-gluon plasma 
(QGP) state might lead to antibaryon abundances greatly exceeding the val- 
ues in a chemically equilibrated hadron gas [5]. However, due to conservation 
of entropy, those abundances may be considerably diluted to a level very close 
to chemical equilibrium [6]. On the other hand, antibaryons have a large an- 
nihilation cross section, which may be modified in the baryon rich colliding 
systems [7]. Measurements of p at the AGS may also contain a large feed-down 
contribution from the decay of antilambdas, A — > p + tt^, as well as from 
other antihyperons (F). Microscopic calculations have given some indication 
of how the annihilation process might change in the dense environment of a 
heavy-ion collision. It has been shown that the amount of annihilation in the 
nuclear medium depends strongly on the formation time of the hadrons [8]. In 
view of all this, it came as a surprise that both in Si+Au and Au-|-Au collisions 
the measured p/p ratio is well described by computing the yields from a fire- 
ball in complete thermal and chemical equilibrium [9,10]. Since the conditions 
for thermal equilibrium are equal to those for hydrodynamic evolution of the 
fireball, it is of great interest to determine experimentally the flow features for 
antiprotons. 

The hydrodynamic evolution of the fireball leads to characteristic flow pat- 
terns such as radial expansion [9] and azimuthal anisotropics, first discovered 
in ultra-relativistic nuclear collisions by ref. [11]. Since then, the study of 
anisotropics in the azimuthal distribution of particles, also called anisotropic 
(directed, elliptic, etc.) transverse flow, plays an important role in high en- 
ergy nuclear collisions [12-14]. Detailed measurements (for surveys see refs. 
[15,16]) have revealed that in a non-central high energy heavy-ion collision, 
nucleons in the backward hemisphere arc preferentially emitted in the direc- 
tion of the impact parameter vector pointing from the target to the projectile. 
Antiprotons co-moving with those nucleons have a greater probability of anni- 
hilation and rescattering, and this could result in an anticorrelation with the 
nucleon directed fiow - the so-called anti-flow of antiprotons in nuclear colli- 
sions predicted by Jahns et al. [17,18]. Measurements of the directed fiow of 
antiprotons are important for understanding the role of annihilation in dense 
nuclear matter. They could also provide insight into the mechanism of anti- 
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quark production in heavy-ion collisions. A surprising result reported in refs. 
[19,20]. is the possible anomalous enhancement of the K/p ratio even beyond 
values expected for complete chemical equilibration [9]. One possible explana- 
tion could be that A baryons have a lower annihilation cross section than p. 
A straightforward way to test this hypothesis would be the measurement of p 
and A (anti-)flow [21]. 

Antiproton production in heavy-ion induced reactions at the AGS has been 
studied experimentally by the E858/E878 [22,23], E814/E877 [24,25], E866 
[26], E886 [27], E864 [28,19], E802 [29,30], and E917 [20] collaborations. In 
the current paper we present the first results on the measurements of directed 
flow of antiprotons detected in the E877 spectrometer in Au+Au collisions at a 
beam momentum of 11. 5A GeV/ c. The results were obtained from the analysis 
of 46 millions of central Au-|-Au events recorded during the 1995 heavy-ion 
run. 



2 E877 apparatus 

Fig. 1 shows the E877 experimental setup which is also discussed in [11,31,32]. 
The trajectory of each beam particle was measured by two silicon micro-strip 
beam vertex detectors (BVER), located at 2.8 m and 5.8 m upstream of the 
target. For the 1995 run, each detector was upgraded from single-sided silicon 
wafers with one-dimensional pitch of 50 /xm to double-sided wafers with a 
200 iim pitch in both the x and y directions [33]. Using these detectors the 
coordinates and angle of beam particles at the target were determined with 
an accuracy of 300 /im in position and 60 //rad in angle. 

The determination of the centrality of the collision and of the reaction plane 
orientation were done using the transverse energy distribution measured in 
the target calorimeter (TCAL), and participant calorimeter (PCAL). Both 
calorimeters had 27r azimuthal coverage and, combined, provided nearly com- 
plete polar angle coverage: TCAL and PCAL covered the pseudorapidity re- 
gions —0.5 < ?7 < 0.8 and 0.8 < 77 < 4.2, respectively [31]. 

Charged particles emitted in the forward direction and passing through a col- 
limator (—134 mrad< Ohorizontai < 16 mrad, —11 mrad< OverUcai < 11 mrad) 
were analyzed by a high resolution magnetic spectrometer with a horizontal 
bend-plane. The spectrometer acceptance covered mostly the forward rapidity 
region. The momentum of each particle was measured using two drift cham- 
bers, DC2 and DC3, whose pattern recognition capability was aided by four 
multi-wire proportional chambers (MWPC). The average momentum reso- 
lution was Ap/p 3%, limited by multiple scattering. A time-of-flight ho- 
doscope (TOFU) located directly behind the tracking chambers provided the 
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time-of-flight information with an average resolution of 85 ps. Energy loss 
measurements in TOFU and in a Forward Scintillator (FSCI) array located 
approximately 30 m downstream of the target were used to determine the 
particle charge. 

For the 1995 run, the spectrometer was upgraded with upstream tracking 
chambers. Two identical multi-wire proportional chambers with highly seg- 
mented chevron-shaped readout pads, VTXA and VTXB, were instrumented 
and placed at 2 m and 2.25 m downstream of the target, just in front of the 
spectrometer magnet. They provided a precise measurement (about 300 /im 
resolution) of the track coordinate in the bending plane of the spectrometer 
before deflection in the magnetic fleld. A description of the design, implemen- 
tation and performance of the vertex detectors can be found in [34]. 



3 Antiproton identification 

The particle identiflcation was performed by combining measurements of mo- 
mentum, velocity, and charge of the particle. To reduce the background in 
the identification of antiprotons, tracks were required to have a confirmation 
from both upstream VTX detectors. Tracks were also required to have com- 
plete information from the cathode pad readout in both drift chambers. Fig. 2 
shows the distribution of particle mass squared, calculated from the measured 
momentum and velocity, for negative particles with momentum less than 4 
GeV/c. A clear antiproton peak is seen, with good signal-to-background ra- 
tio. The contribution of the background increases with momentum (see lower 
panel of Fig. 2) and becomes dominant at p > 4.5 GeV/c. A maximum mo- 
mentum of 4 GeV/c was required for clean identiflcation of antiprotons. The 
mass resolution is quantitatively understood in terms of the intrinsic detector 
resolution, multiple scattering and time-of-flight resolution. For a positive an- 
tiproton identiflcation we selected, in a plot of momentum versus the measured 
mass, a region of ±2.0o"m2(p) around the antiproton mass peak, where (T„i2(p) 
represents the resolution in the mass squared at a given particle momentum 
p. The dependence of the mass resolution on momentum was assumed to be 
equal to that of protons. 

Fig. 3 shows the antiproton acceptance in the transverse momentum pt and 
rapidity y coordinates. The E877 spectrometer mainly covers the low pt region 
{pt < 0.4 GeV/c) forward of ycm = 1-6 (1.6 < y < 2.2). The flnal data sample 
comprises about 750 antiprotons with about 400 in the 1.8 < y < 2.2 rapidity 
region. 

Since tracking only starts 2 m downstream of the target, the E877 spectrom- 
eter is not well suited for separation of primary antiprotons from those fed 
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down from anti-hyperon decays. Because of the large mass asymmetry and 
small decay momentum in these decays, the antiproton track points back to 
the target. A previous extensive study [35] based on a Monte-Carlo simulation 
showed that antiprotons are reconstructed for 70% of the A and for 100% of 
the E~ decaying such that the p is emitted into the acceptance of the spec- 
trometer and consistent with originating from the target. 



4 Results and Discussion 



4-1 Experimental data 



The azimuthal anisotropy in particle production is studied by means of Fourier 
analysis of azimuthal distributions [11,31,32,36,37] with respect to the reac- 
tion plane. We study the rapidity, transverse momentum, and centrality de- 
pendence of the Fourier coefficients v!^ (amplitude of n-th harmonic) in the 
decomposition: 

(PN 1 d'^N 

E—— = — -(1 + cos + 2^;^ cos 20 + ...), 

d-^p 27r ptdptdy 



where the azimuthal angle = (piai, — ipr is taken with respect to the reaction 
plane orientation. Directed flow is quantified by the dipole Fourier coefficient 
v[. 

Similarly to the analysis performed in [31,32], the reaction plane angle is 
determined from the measured transverse energy distribution in four non- 
overlapping pseudorapidity windows. The reaction plane resolution, i.e. the 
accuracy with which the reaction plane orientation is determined, is evalu- 
ated by studying the correlation between flow angles measured in different 
windows. Finally, the flow signals are corrected for the flnite reaction plane 
resolution. Details of this procedure are described in [31,32]. In the following, 
only coefficients vi corrected for the reaction plane resolution will be shown. 

The experimental acceptance for p is close to mid-rapidity (Fig. 3), where 
the directed flow changes sign. To optimize a possible flow signal, we skip 
the region very close to mid-rapidity and present, in Fig. 4, the measured 
azimuthal distributions of antiprotons emitted at rapidity y > 1.8, for two 
collision centralities. To determine Vi, the distributions were fitted with a dis- 
tribution /(0) = l + 2vi cos (j). A pronounced minimum is observed at = for 
semicentral collisions {a = 10 — 26% a^eo), indicating that the antiproton pro- 
duction is strongly anti-correlated with the flow of nucleons. For more central 
collisions (top 10% of ageo), the antiproton azimuthal distribution is, within 
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errors, isotropic {vi ~ 0). No or very little proton flow has been observed for 
the most central Au+Au collisions [32]. 

In order to check the contribution of the background to the observed flow- 
signal we evaluated the azimuthal distributions for the particles whose mea- 
sured mass was outside the window for positive antiproton identification. For 
the background sample we selected, in a plot of momentum versus the mea- 
sured mass, a region of ±2-6 a^2 (p) away from the p mass peak. On the low 
mass side we cut the region less than 3.5 cr^2(p) away from the K~ peak. 
The azimuthal asymmetry for the background sample was found to be con- 
sistent with zero (f i = —0.023 ± 0.053 for 10-26% cTgeo centrality window and 
Vi = 0.021 ± 0.038 for 0-10% (Xgeo centrality window). We also checked the 
t'i(pt)-dependence for the background sample and found no signal within the 
statistical errors for all pt bins. The data shown do not include the correction 
due to background subtraction. Assuming a typical signal-to-background ratio 
of 1.2 and that background has no asymmetry, the average vi for 10-26% ageo 
centrality window increases to ~ — 0.19±0.07. 

The measured transverse momentum dependence of vi is shown in Fig. 5. 
Large negative values of vi are observed for pt > 0.1 GeV/c. It should be 

noted that the absolute value of the proton flow signal at similar rapidities 
and transverse momenta is significantly smaller [32]. The right hand panel 
shows a result consistent with zero flow for the more central bin. 



4.2 RQMD predictions 

The Relativistic Quantum Molecular Dynamics (RQMD) model [38] has been 
widely used in describing relativistic heavy-ion collisions. It combines classical 
propagation of all hadrons with string and resonance excitations in the primary 
collisions of nucleons from target and projectile. Overlapping color strings may 
fuse into so-called ropes [39]. Subsequently, the fragmentation products from 
rope, string and resonance decays interact with each other and with the origi- 
nal nucleons, mostly via binary collisions. These interactions drive the system 
towards equihbration [40] and are responsible for the development of collective 
flow, even in the pre-equilibrium stage. If baryons are surrounded by other 
baryons they acquire effective masses. These effective masses are generated 
by introducing Lorentz-invariant Skyrme-type quasi-potentials into the mass- 
shell constraints for the momenta, which simulates the effect of a mean field 
[41]. There are no potential-type interactions in the so-called cascade mode of 
RQMD. 

For antiprotons, RQMD combines a large enhancement of initial production 
via hadronic multi-step processes with strong absorption through the free p 
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absorption cross section [42], which results in the absorption of a large frac- 
tion of the produced antiprotons (~ 90% of p are reabsorbed for minimum bias 
Au+Au collisions [18]). The newest version of the model (RQMD v2.3) uses 
the formation of a quasi-bound state (proton-antiproton molecule) to signifi- 
cantly reduce the possibility of antiproton annihilation. However, the presence 
of a mean color potential could possibly have a nmch more significant effect 
on the final antiproton distribution [43] . A systematic study of the antibaryon 
production with various projectile-target combinations and at different en- 
ergies would be useful for refining this approach and for understanding the 
annihilation process in relativistic heavy-ion collisions. 

Predictions of the RQMD model run in cascade mode for the prdependence 
of Vi for antiprotons are compared to the data in Fig. 5. For the model, we 
only used antiprotons produced in the experimental acceptance window. The 
magnitude of the negative signal observed in the data for semicentral collisions 
is comparable, within the errors, to that predicted by the model although, 
systematically, the model values are somewhat smaller. 

A current question of interest and debate is the relative importance of in- 
cluding nuclear mean field effects in intra-nuclear cascade models. At lower 
energies, experimental data and a detailed comparison with theoretical pre- 
dictions have been discussed in [14]. An attractive p potential of the order of 
100 to 150 MeV at normal nuclear density is needed to reproduce the size and 
shape of the experimental spectra [44]. It is expected that collective flow may 
be an important observable for evaluating the importance of mean field effects 
[43]. It was shown, in particular that, in semicentral Au+Au collisions, protons 
exhibit a strong azimuthal anisotropy, with the amplitude of flow well repro- 
duced by RQMD only if the effects of m,ean field are included [32]. It should 
be noted, however that, even including the mean fleld, only the p^- integrated 
signal is reproduced by the model, while the predicted p^-dependence signifi- 
cantly differs from the data. In Fig. 6 we present the RQMD v2.3 calculations 
for proton and antiproton mean values of the transverse momentum projected 
onto the reaction plane, {px)-, as a function of rapidity for Au+Au collisions at 
11. 5 A GeV/c with impact parameters 6 < 10 fm. Calculations are presented 
for a pure cascade and mean-field modes of RQMD. Since the calculations are 
very CPU intensive we have a limited antiproton statistics for the mean field 
calculations: about 500 p's compared to more than 8700 p's for the cascade 
mode. The flow of nucleons is enhanced due to additional repulsion in the 
mean fleld approach. For antiprotons, on the other hand, the effective mean 
field potential in nuclear matter is expected to be attractive. The antinucle- 
ons are pulled towards the nucleon source, therefore the anti-flow should be 
considerably weakened [45], as observed in Fig. 6. Such a weakening is also 
supported by the results of Spieles et al. [43], where it was shown that, when 
the infiuence of the real part of an antinucleon-nucleus optical potential on 
the p momentum spectra is included into the calculations, the antiproton an- 
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nihilation is indeed suppressed in the heavy- ion environment, such that the 
sum of in-medium annihilation cross section and clastic cross section is closer 
to the free annihilation cross sections which are employed in the cascade mode 
of RQMD v2.3. On the other hand, the amplitude of the measured antiproton 
flow signal is rather well described by the pure cascade calculations (see Fig. 5). 
A similar observation has been recently made by comparing the measurements 
of directed flow with the cascade and mean fleld RQMD predictions [46], 
where it has been shown that the cascade calculations reproduce the data 
better than the mean field calculations. Further measurements and detailed 
calculations in the mean field mode are required and this question cannot be 
fully settled until antiproton yield and fiow have been measured over a larger 
phase space. 



5 Conclusion 

We have observed, for the first time, a strong azimuthal anisotropy in an- 
tiproton production in semicentral Au-|-Au collisions which is anti-correlated 
with the nucleon emission. The observed large anti-flow for p indicates that 
a strong annihilation process is involved in the dense nuclear matter as pre- 
dicted by previous theoretical studies. The presented data have confirmed the 
need for a thorough investigation, both experimentally and theoretically, of 
antibaryon production in nucleus-nucleus interactions. Whereas the measured 
antiproton yields/spectra can be described within a thermal model assum- 
ing a chemical/thermal equilibration and collective expansion of the system, 
antiprotons exhibit a directed flow which is even stronger than expected by 
the models incorporating mean-field effects. A complementary measurement 
of the directed fiow of antihyperons such as A in the future combined with 
our results would provide necessary information for the decoupling of the fiow 
contribution of primary antiprotons and that of antihyperons. Furthermore, 
the comparison between the directed flow behavior of these two contributions 
could give crucial information for understanding the observed unexpected A/p 
ratios. 
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Fig. 1. E877 experimental setup. 
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Fig. 2. Upper plot: the mass squared distribution in the antiproton peak region 
for negative particles with momenta less than 4 GeV/c. The Lower plot shows the 
same distribution for three different momentum intervals. The solid line is a fit by 
a Gaussian and an exponential. 
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Fig. 3. Antiproton acceptance in the {pt,y) phase space covered by the E877 spec- 
trometer. Midrapidity is y = 1.6. 
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Fig. 4. Antiproton azimuthal distributions measured for rapidity 1.8 < y < 2.2 and 
two collision centralities. The distributions are normalized such that the average is 
unity. The solid line is a fit by a 1 + 2vi cos (f) distribution. 
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Fig. 5. vi{pt) of antiprotons for rapidities 1.8 < y < 2.2 and two collision ccntralitics. 
The data (solid circles) are compared with the calculations using the RQMD v2.3 - 
cascade mode (open circles). 
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Fig. 6. Mean values of the transverse momentum projected into the reaction plane 
ipx) of protons (dashed lines) and antiprotons (solid lines) as a function of rapidity 
for Au+Au collisions at 11. 5 A GeV/c with impact parameters 6 < 10 fm generated 
by using the RQMD v2.3 {cascade mode - thin lines; mean field mode - thick lines). 
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